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ABSTRACT 

We present a study of the 2014 January 6 solar energetic particle event which produced a small ground level 
enhancement (GLE), making it the second GLE of this unusual solar cycle 24. This event was primarily observed 
by the South Pole neutron monitors (increase of ~2.5%) while a few other neutron monitors recorded smaller 
increases. The associated coronal mass ejection (CME) originated behind the western limb and had a speed of 
1960 km s -1 . The height of the CME at the start of the associated metric type II radio burst, which indicates the 
formation of a strong shock, was measured to be 1.61 Rs using a direct image from STEREO -A /EUVI. The CME 
height at the time of the GLE particle release (determined using the South Pole neutron monitor data) was directly 
measured as 2.96 Rs based on STEREO-A /CORl white-light observations. These CME heights are consistent 
with those obtained for GLE7 1 , the only other GLE of the current cycle, as well as cycle-23 GLEs derived using 
back-extrapolation. GLE72 is of special interest because it is one of only two GLEs of cycle 24, one of two 
behind-the-limb GLEs, and one of the two smallest GLEs of cycles 23 and 24. 
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1. INTRODUCTION 

Solar energetic particle (SEP) events, where particles accel- 
erated to GeV energies are subsequently detected on the ground 
as a result of the air-shower process, are known as ground level 
enhancements (GLEs). With a typical detection rate of a dozen 
GLEs per cycle, an average of 16.3% SEP events were GLEs 
in cycles 19-23 (Oliver et al. 1982; Oliver 2006; Shea & Smart 
2008; Mewaldt et al. 2012; Nitta et al. 2012; Gopalswamy et al. 
2012a). In cycle 24, this fraction is much smaller (6.4%) with 
2 GLEs out of 31 large SEP events (Gopalswamy et al. 2014). 
This is also much smaller than the ratio of 18% obtained when 
the first five years of cycle 23 are considered. GLEs are typi- 
cally associated with intense flares (median soft X-ray intensity 
~X3.8) and fast coronal mass ejections (CMEs; average CME 
speed ~2000 km s _1 ; see Gopalswamy et al. 2012a). 

In cycle 24, only one GLE (2012 May 17; see Gopalswamy 
et al. 2013b) occurred prior to the 2014 January 6 GLE (GLE72) 
studied here. The 2014 January 6 GLE is unique because despite 
being a backside event, observations from the Solar Terrestrial 
Relations Observatory (STEREO) provide a view of the source 
region and the early CME propagation in full detail. This is the 
first set of observations that covers the spatial domain of interest 
for a behind-the-limb GLE. We exploit this unprecedented 
opportunity to understand the early kinematics of the associated 
CME and shock formation based on extreme ultraviolet (EUV) 
measurements and the height of the CME at the time of particle 
release. 

In GLEs, the first arriving particles can be considered to 
be relatively scatter-free, and so they provide information on 
the solar particle release (SPR) time, i.e., the time when the 
accelerated particles were injected into interplanetary space. It 
has been widely accepted that the majority of GLE particles are 
accelerated by CME-driven shocks (Reames 2013; Reames et al. 
2013). Thus, the knowledge of the SPR time together with the 
coronagraphic observations can be used to determine the height 
at which the GLE particles are released. For a shock driven by a 
high-speed CME, the height of the shock can be approximated 


as the height of the CME itself. Combined measurements from 
the Large Angle and Spectrometric Coronagraph (LASCO) 
on board the Solar and Heliospheric Observatory (SOHO), 
with a field of view (FOV) of 2.5-32 Rs (Brueckner et al. 
1995) and the inner coronagraph (CORl) on STEREO with 
a FOV extending closer to the solar surface (1.4 Rs; Thompson 
et al. 2003), provide accurate CME kinematics, the height of 
shock formation, and the height of the shock at the SPR time. 
Additional details of the eruption, down to the solar surface, 
are obtained from the Extreme Ultraviolet Imager (EUVI) on 
STEREO (Wuelser et al. 2004; Howard et al. 2008). It is 
important to understand the GLE of 2014 January 6, not only 
because it is the only other GLE of the unusually weak cycle 
24 but also because there have been nearly 40 large (flare soft 
X-ray magnitude ^ M5) SEP events in this cycle, 8 of which 
were well-connected (source region longitude range: W55- 
W88) but did not produce a GLE (Gopalswamy et al. 2013a). 

2. OBSERVATIONS AND DISCUSSION 

As shown in Figure 1(a), the South Pole neutron monitor 
(SPNM; Bieber et al. 2013; Evenson et al. 2011) recorded the 
onset of GLE72 around 07:58 UT and measured a maximum 
rise of ~2.5% above the cosmic-ray background. A few other 
neutron monitors, including McMurdo and Tixie Bay, observed 
smaller increases (<2%), making it difficult to determine the 
anisotropy during GLE72. This event was observed by parti- 
cle detectors on board the Geostationary Operational Environ- 
mental Satellites (GOES) at higher energies (>700 MeV), and 
by detectors on the Advanced Composition Explorer (ACE), 
STEREO, and Wind spacecraft (von Rosenvinge et al. 1995) at 
lower energies (^100 MeV). The GOES 420-510 MeV proton 
channel observed the onset at ~08: 12 UT. The increase in inten- 
sity of >700 MeV protons from GOES (Figure 1(b)) confirms 
the 2014 January 6 event to be a GLE because all GLEs in 
cycles 23 and 24 have discernible enhancements in the GOES 
>700 MeV energy channel (Gopalswamy et al. 2014; Tylka & 
Dietrich 2009). Both GOES-13 and GOES-15 (Rodriguez et al. 
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Figure 1 . Left: (a) count rates from the South Pole neutron monitor (http://neutronm.bartol.udel.edu). A rise of ~2.5% was observed, (b) Proton intensity time profile 
from GOES/ EPEAD. GOES observed a ~67% increase in the >700 MeV channel. The onset at SPNM (dotted vertical line) occurs before GOES. The normalized 
solar particle release time was 07:55 UT. Right: (upper) m-Type II radio burst observed by the Learmonth radio observatory (www.ips.gov.au/Solar/3/1) starting at 
07:45 UT (the arrows point to the fundamental-harmonic (F-H) bands) and (lower) the decameter-hectometric (DH; pointed by an arrow) to kilometric type II observed 
by Wind /WAVES. 


2014) detected this event, but the GOES- 13 proton fluxes were 
greater at higher energies. 

The Learmonth radio observatory observed a metric type II 
(m-Type II) radio burst at 07:45 UT, which lasted until 08:05 
UT (Figure 1, upper right). The radio and plasma wave receiver, 
WAVES (Bougeret et al. 1995), on the Wind spacecraft detected 
the type II burst in decameter-hectometric wavelengths (DH- 
Type II) from 14 MHz down to 1 MHz and in the kilometric 
domain (km-Type II) from 1 MHz to 200 kHz (Figure 1, lower 
right). A type II burst extending from metric to kilometric 
wavelengths is an indication that a CME is able to drive a 
strong shock from the low solar corona into the interplanetary 
medium and is typical of events with GLEs (Gopalswamy 2011; 
Gopalswamy et al. 2010). 

The CME of interest was observed by STEREO- A/ COR1, 
STEREO- A /EUVI, and SOHO/ LASCO (Figure 2). STEREO-A 
(located at W151) confirmed that the CME originated behind 
the limb (Figure 2(a)). STEREO-A/ EUVI, STEREO -A /COR1, 
and SOHO /LASCO-C2 first observed the CME at 07:30 UT, 
07:55 UT, and 08:00 UT, respectively. These multi-spacecraft 
coronagraphic observations provide excellent data to character- 
ize the CME kinematics. The eruption was also observed by the 
Atmospheric Imaging Assembly (AIA) on the Solar Dynamics 
Observatory ( SDO\ Pesnell et al. 2012; Lemen et al. 2012) over 
the west limb. A C2.2 solar flare starting at 07:30 UT was re- 
ported by the online Solar Geophysical Data, but this flare class 
does not represent the true magnitude of the soft X-ray flare 
because most of the flare structure was occulted by the west 
limb. 

3. ONSET OF THE EVENT AND 
CALCULATION OF SPR TIME 

We used the two-minute data from SPNM to estimate the 
SPR time. Since this is a small GLE with some fluctuations 
in the background intensity just before the onset of the event, 
we estimate the onset of the GLE to lie between 07:54 UT 
and 08:02 UT. We take the median of this, 07:58 UT, as the 


onset time for GLE72 at SPNM. The solar wind speed was 
~ 350-400 km s _1 . The geomagnetic conditions were quiet with 
a Kp index of +1 around the event onset (Geomagnetic data 
service, Kyoto: http://wdc.kugi.kyoto-u.ac.jp). The asymptotic 
directions for several neutron monitors are shown in Figure 3. 
The direction of the interplanetary magnetic field (IMF) fluctu- 
ated near the event onset. The particle arrival directions based 
on the IMF are also shown in Figure 3. We assume that after be- 
ing accelerated and injected into the interplanetary medium, the 
particles traveled a nominal path length of 1 .2 AU before reach- 
ing Earth. SPNM has a geomagnetic rigidity cutoff of 0.09 GV 
(~0.004 GeV), but atmospheric effects dominate the NM mea- 
surements. The high altitude of SPNM (2820 m) can allow 
observation of lower energy particles (~300 MeV). However, 
the first arriving GLE particles would be of the highest energies 
from the accelerated population. In addition, the fact that GOES 
recorded an onset of >700 MeV protons several minutes after 
the event onset at SPNM also confirms that the GLE72 onset 
at SPNM was due to particles of energies >700 MeV. Thus, 
we take 1 .0 GeV to be the energy of the first arriving particles 
at SPNM. The 1.0 GeV protons (speed of ~2.63x 10 5 km s -1 ) 
would take 11.4 minutes to travel a distance of 1.2 AU, and 
therefore the SPR time is estimated to be 07:47 UT. In the 
CME-normalized frame of reference (incorporating ~8.3 min- 
utes for the propagation of solar electromagnetic radiation to 
Earth) the SPR time is ~07:55 UT. This timing is consistent 
with many previous observations which show that the shock 
forms several minutes before SPR occurs (Kahler et al. 2003; 
Gopalswamy et al. 2012a, 2013a). 

4. CME KINEMATICS AND HEIGHT AT SPR 

We utilized observations from STEREO and SOHO/ LASCO 
to obtain the early kinematics of the CME using the graduated 
cylinder model of a flux-rope (Themisien 2011). The flux-rope 
in this model has conical legs with a circular cross section and a 
circular front, and it expands self- similarly. The flux-rope fit to 
the CME observations provides the height-time information and 
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Figure 2. Left: (a) 2014 January 6 CME observed by STEREO-A / EUVI at the onset of m-Type II (07:45 UT) burst, (b) The post-flare loops seen above the west limb 
in the SDO / AIA 193 A image (07:48 UT). (c) STEREO-A/ COR1 and EUVI image of the CME at SPR (07:55 UT). (d) The CME as seen by LASCO/C2 (08:12 UT). 
Right: CME kinematics obtained by fitting a flux-rope to STEREO-A/ EUVI, COR1, COR2, and SOHO / LASCO observations. The height and acceleration of the 
CME were ~2.96 Rs and 2.08 km s -2 , respectively, at SPR. The CME speed at SPR was 1658 km s _1 . The GOES 1. 0-8.0 A soft X-ray flux is superposed. 


Asymptotic directions for neutron monitors during GLE72 
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CS: Cape Schmidt MA: Mawson NO: Norilsk TH: Thule 

Figure 3. Neutron monitor viewing directions (red) at 07:58 UT on 2014 January 6. The particle arrival directions based on the IMF direction are also shown (blue). 


width for the CME as well as the source region coordinates. The 
height- time, speed, and acceleration for the GLE72 CME are 
depicted in Figure 2. The fit determined W 102 to be the source 
longitude of the flux-rope. The CME reached its peak speed 
(~1960 s -1 ) and experienced peak acceleration (~2.08 km s -2 ) 
at 08:00 UT and 07:55 UT, respectively, and thus the CME 
acceleration was at its maximum at the time of SPR. The 
drag force decelerated the CME and the acceleration ceased 
between 08:00 and 08:05 UT. The time of shock formation 
is indicated by the onset of the m-Type II burst (Reames 
2013; Gopalswamy et al. 2013b, 2012a, 2012b; Reames 2009a, 
2009b). Type II onset times provide crucial information about 


particle acceleration because electromagnetic radiation does not 
depend on magnetic connectivity. STEREO-A / EUVI observed 
the associated eruption at 07:45 UT (Figure 2(a)) during the 
onset of m-Type II, and hence provides a direct measurement of 
the CME height at the shock formation as 1.61 Rs. STEREO-A/ 
COR1 observed the CME at 07:55 UT (Figure 2(c)), thus giving 
a CME height of 2.96 Rs at SPR. 

The timeline of events during the GLE based on the mea- 
surements and calculations is shown in Table 1 . The CME ac- 
celerated from 07:30 UT. After its formation at 07:45 UT, the 
CME-driven shock traveled for ~ 10 minutes (acceleration time) 
before particles were released into the interplanetary medium. 
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Table 1 

Timeline of Events for GLE72 


Event Time 

(UT). 

CME onset 07:30 

CME first observation by STEREO-A / EUVI 07:30 

Metric type II onset 07:45 

Normalized SPR time 07:55 

CME first observation by STEREO-A /CORA 07:55 

GLE particles arrival at SPNM 07:58 

CME first observation by SOHO / LASCO 08:00 


Note. All times are Earth times and propagation times for electro- 
magnetic emissions have been considered in this table as explained 
in the text. 

5. COMPARISON OF GLE72 WITH OTHER GLES 

We now compare the 2014 January 6 GLE with a few GLEs 
in cycles 23 and 24 which have some overlapping properties. In 
particular, we consider the GLEs with longitude and intensity 
similar to those of GLE72. 

5.1. Comparison with the Solar Cycle 24 GLE (GLE71) 

GLE71, almost a magnitude higher than GLE72, was as- 
sociated with an M5.1 flare from W89 and a fast CME 
(~1997 km s _1 ). The CME was observed by STEREO-A/ 
COR1, giving the CME height at SPR to be 2.32 Rs. The ac- 
celeration of the CME at the time of shock formation and SPR 
were 1.77 km s -2 and 1.51 km s -2 , respectively (Gopalswamy 
et al. 2013a). We note that the CME speeds were comparable 
for GLE71 and GLE72 but the CME acceleration at SPR was 
slightly higher for GLE72. In addition, the GOES proton time 
profiles in Figure 1 are comparable to those in Figure 3(d) of 
Gopalswamy et al. (2013a), with similar intensity enhancements 
in the >700 MeV channel (~84% for GLE71 versus ~67% for 
GLE72). 

5.2. Comparison with Solar Cycle 23 GLEs 

Here, we compare the GLE72 CME heights at the m-Type II 
onset and SPR, and the acceleration at SPR with those of the 
cycle 23 GLEs. Figure 4 shows the CME heights at the time 
of shock formation and SPR for GLEs of solar cycle 23 as a 
function of the source longitude. The CME height at SPR ( h s 
in Rs) shows a parabolic dependence on the source longitude 
(A in degrees): h s = 2.55 + [(A.-51)/35] 2 (Gopalswamy et al. 
2012b). Reames (2009a, 2009b) obtained a similar figure using 
an independent technique. GLE72 lies ~37% below this curve. 
We note that the source longitude distance from the nominal 
magnetic footpoint of Earth were similar for GLE 65 (2003 
October 28, longitude: E02) and GLE72. The GLE65 flare 
magnitude was XI 7, the most intense of all GLEs of cycles 
23 and 24. The CME speed and acceleration for GLE65 were 
2754 km s -1 and 1.53 km s -2 , respectively, and the shock formed 
at a relatively low height of 1.27 Rs. The CME height at SPR 
is closer to the parabola for GLE65. The CME heights at SPR 
exhibit a standard deviation of ~1.61 Rs from the parabola. 
The GLE72 CME is 1.71 Rs below the parabola but a few GLE 
CMEs deviate farther (e.g., GLE62 CME is ~4.61 Rs higher). 

In an extensive study of the CMEs associated with cycle 
23 GLEs, Gopalswamy et al. (2012a) listed the mean and 
median CME accelerations to be 2.3 km s -2 and 2.02 km s -2 , 
respectively. We note that the 2014 January 6 GLE CME has 
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GLE and m-Type II Height vs. Longitude 



Source Longitude [deg] 


Figure 4. Height of the CME at shock formation and SPR plotted against the 
source longitude of cycle 23 GLEs. The CME heights at SPR can be fitted with 
a parabola as shown. The 2014 January 6 event (GLE72) as well as GLEs61, 
65, 68, and 71 are shown superposed. The green and blue diamonds indicate 
heights determined from linear and quadratic extrapolations. 

an acceleration of 2.08 km s -2 at SPR, and thus agrees with the 
acceleration of CMEs associated with cycle 23 GLEs. 

5.3. Comparison with a Small GLE of Solar Cycle 23 (GLE68) 

Because the 2014 January 6 GLE was a small event with a 
maximum increase of ~2.5%, we compared it to the smallest 
GLE of cycle 23, that is, the GLE of 2005 January 17 (GLE68), 
for which a maximum increase of ~3% was observed. As 
shown superposed on Figure 4 (GLE68 marked in black), the 
CME heights at shock formation and SPR time were 1.44 Rs 
and 2.72 Rs, respectively. The CME had a maximum speed of 
2802 km s -1 and an acceleration of 2.92 km s -2 (Gopalswamy 
et al. 2012a). These are higher than those for GLE72 but the 
CME heights at shock formation and SPR are still consistent. 
It is interesting that although the GLE68 source location, W25, 
was within the well-connected solar region and the eruption was 
associated with a fast CME with a large acceleration, it produced 
a very small GLE. One possible reasons for this is the poor 
latitudinal connectivity for GLE68 (Gopalswamy et al. 2013a). 
Even though the flare location was N14W25, the flux-rope fit 
indicated a location of N18W25. The latitudinal connectivity 
became worse because of the unfavorable solar B0 angle (~4?7), 
making the effective location N23W25. The latitudinal distance 
to the ecliptic (23°) is 10° greater than the nominal latitudinal 
distance (13°) for cycle-23 GLEs. The CME height at SPR falls 
on the parabola for GLE68. 

5.4. Comparison with a Behind- the -limb GLE (GLE61 ) 

The only backside GLE of cycle 23 is the 2001 April 18 GLE 
(GLE61). The solar source was at W117. The associated CME 
had a speed of 2712 km s -1 and an acceleration of 2.38 km s -2 . 
In this event, the time delay between the onset of m-Type II and 
SPR was about 15 minutes (Gopalswamy et al. 2012a). The time 
delay for GLE72 is about 10 minutes. As shown in Figure 4, the 
GLE61 CME was at heights of 1.47 Rs and 6.01 Rs at the 
m-Type II onset and the SPR time, respectively. This is 
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consistent with a particle event originating behind the limb, be- 
cause the particle release occurs at a later time when the shock 
crosses the magnetically well-connected field lines. These two 
GLEs are on the right arm of the parabola because both oc- 
curred behind the west limb. GLE61 was ~27° behind the limb 
compared to only ~12° behind the limb for GLE72. Therefore, 
the CME height at SPR is smaller for GLE72, as expected. 
It is not clear whether or not the smaller height at SPR for 
GLE72 has anything to do with the state of the heliosphere 
during a weak cycle. The 2014 January 6 GLE and the 2001 
April 1 8 GLE form a very interesting subset of GLEs that pose 
the greatest challenge to the flare acceleration mechanism for 
GLE particles. In the case of acceleration by flare, the parti- 
cles would have been accelerated by a narrow source at the 
reconnection region, typically at a height of ^0.05 Rs above 
the solar surface. Therefore, a flare acceleration region would 
not only be occulted, but would also be poorly connected to 
Earth (Gopalswamy et al. 2012a), and thus would not result in 
a GLE. However, in both of these GLEs (61 and 72), the ener- 
getic particles were observed promptly at Earth. In the scenario 
of shock acceleration, the CME-driven shock expands and the 
GLE particles can be observed when the eastern flank crosses 
the field lines that are magnetically connected to Earth. There- 
fore, flare accelerated particles would need to interact with the 
CME flux rope to reach the well-connected field lines. In this 
case, the first arriving particles will not remain scatter-free, 
contradicting the observed high anisotropy at the beginning of 
the GLEs. 

6. SUMMARY 

We have analyzed the GLE of 2014 January 6 using STEREO - 
A/ EUVI, STEREO- A/ COR1, and SOHO/ LASCO measure- 
ments to obtain accurate CME kinematics. GLE72 is unique 
in many aspects as it falls into different categories of a few spe- 
cial GLEs in cycles 23 and 24. In the unusual cycle 24, where 
there has been a dearth of GLEs, this is one of the only two 
GLEs. This is also one of the two smallest GLEs (GLE68 and 
GLE72) of cycles 23 and 24. In addition, GLE72 is one of the 
only two behind-the-limb GLEs (GLE61 being the other one) 
of cycles 23 and 24, thus posing a challenge to flare acceleration 
theory for GLE particles. Our analysis of GLE72 is summarized 
as follows. 

1. The observation of >700 MeV protons by GOES confirms 
the 2014 January 6 SEP event as a GLE, for which a small 
intensity enhancement was observed by the South Pole and 
a few other neutron monitors. This is consistent with the 
findings by Gopalswamy et al. (2014) that all of the GLEs 
of cycle ^ 23 and 24 had discernible enhancement in the 
GOES >700 MeV proton channel. 

2. The CME height at the time of shock formation (onset time 
of m-Type II radio emission) was measured directly by 
STEREO-A/ EUVI to be 1.61 Rs. The CME height at SPR 
was measured directly by STEREO-A /CORA to be 2.96 
Rs. Our analysis is consistent with particle acceleration 
by a CME-driven shock. This is the first behind-the-limb 
GLE for which STEREO provided direct CME height 
measurements both at the shock formation and SPR. This 
event compares well with the 2012 May 17 GLE in terms 
of GOES observations, CME heights at shock formation, 
and SPR. 


3 . The CME height at the time of shock formation is consistent 
with that of other GLE- associated CMEs of cycle 23. The 
height of the CME at the SPR time varies with the longitude 
of the source of the GLEs such that the CME height at SPR 
increases as the source longitude goes beyond the best- 
connected regions (W20-80). The CME height at SPR for 
GLE72 is about 37% below the parabola fitted to the cycle 
23 GLE CMEs. 

4. The CME heights at the shock formation and SPR for the 
small GLE72 are comparable with those for the smallest 
GLE of cycle 23 (GLE68). 

5. The CME height at SPR for GLE72 is consistent with that 
for GLE61; these GLEs originated 12° and 27° behind the 
limb, respectively. 
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